PRECURSORS OF UCHII REGIONS & 

THE EVOLUTION OF MASSIVE OUTFLOWS 



Henrik Beuther 

Harvard-Smithsonian Center for Astrophysics 
60 Garden Street 
Cambridge, MA 02138 
USA 

hbeuther@cfa.harvard.edu 



Debra Shepherd 

National Radio Astronomy Observatory 

Socorro, NM 87801 

USA 

dshepher@nrao.edu 



Abstract Since this contributions was meant to cover two subjects which are both in the 
field of massive star formation but which in its details can be discussed sepa- 
rately, this paper is divided in two sections. First, we present characteristics of 
precursors of UCHll regions and their likely evolutionary properties. The second 
section discusses massive molecular outflows, their implications for high-mass 
star formation, and a possible evolutionary sequence for massive outflows. 
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1. Precursors of UCHu regions 

Evolutionary scenarios suggest that High-Mass Starless Cores (HMSCs) 
represent the earliest evolutionary stage of massive star formation (e.g., Evans 
et al., 2002). The next observable evolutionary stage is when a High-Mass Pro- 
tostellar Object (HMPO) forms in the core producing strong millimeter contin- 
uum and mid-infrared emission but no detectable centimeter continuum due to 
free-free emission from ionized gas (e.g., Molinari et al., 1996, Sridharan et al., 
2002). The lack of detectable free-free emission may be because the protostar 
has not reached the main sequence or because the accretion rate onto the stel- 
lar surface is high enough to quench a developing HII region (e.g., Walmsley, 
1995, Churchwell, 2002). Hot Cores are considered to belong to the HMPO 
stage. Soon after the hot core is formed, the central massive protostar produces 
adequate ionizing radiation to form a hypercompact, unresolved, most likely 
optically thick Hn region (HCHu). The HCHu region may still be partially 
quenched by infalling gas which can hinder the expansion of the hypercompact 
Hil region (e.g., Walmsley, 1995, Keto, 2003). Eventually, the hypercompact 
Hil region begins to expand forming the well-studied ultracompact Hn region 
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(e.g., Wood and Churchwell, 1989, Kurtz et al., 1994) and later a more evolved 
Hll region. Theoretical model calculations estimate an evolutionary time-scale 
for massive star formation of about 10 5 yrs, assuming that massive stars form 
via accretion-based processes similar to low-mass star formation (McKee and 
Tan, 2002). A key difference between low- and high-mass star formation is that 
massive stars reach the main sequence while still actively accreting; low-mass 
stars end their accretion phase before they reach the main sequence. Here, we 
discuss the properties of the precursors of ultracompact Hn regions, namely 
HMPOs and HMSCs. 

1.1 High-Mass Protostellar Objects (HMPOs) 

There are several review articles that discuss the earliest stages of massive 
star formation and cover different aspects of the topic in more depth than pos- 
sible here (e.g., Garay and Lizano, 1999, Kylafis and Pavlakis, 1999, Kurtz 
et al, 2000, Churchwell, 2002, Cesaroni, 2004). Therefore, we focus on the 
basic characteristics of HMPOs and discuss a few recent results in this field. 

Basic characteristics. The far-infrared colors of HMPOs are similar to 
those of UCHu regions (e.g., Wood and Churchwell, 1989, Sridharan et al, 
2002). This is mainly due to the fact that they resemble each other with re- 
gard to their gas/dust temperatures and the presence of dense gas. A major 
observational difference is that HMPOs show no or only weak cm continuum 
emission, thus the central objects have not yet produced a significant ionized 
region to trigger the necessary free-free emission. While most of the UCHu 
region luminosity is due to a central H-burning star, a large fraction of the 
HMPO luminosity is still expected to be due to accretion (e.g., Sridharan et al., 
2002). Furthermore, the average observed NH3(1,1) linewidth toward HM- 
POs is Av ~ 2.1km/s compared to average values toward UCHu regions of 
At; ~ 3.1 (Sridharan et al., 2002). This difference implies less turbulent mo- 
tions at the HMPO stage. 

Furthermore, often H2O and/or Class 11 CH3OH maser emission is observed 
toward HMPOs (e.g., Walsh et al, 1998, Beuther et al., 2002d, Codella et al., 
2004). The existence of one or the other type of maser emission is widely re- 
garded as a good signpost for massive star formation (in spite of H2O maser 
emission being also observed toward low-mass star-forming regions), but nev- 
ertheless, there is no general agreement about the physical processes generat- 
ing the maser emission. There is indicative evidence that both maser types are 
produced either within accretion disks or molecular outflows (e.g., Norris et al., 
1998, Walsh et al., 1998, Torrelles et al, 1998, Minier et al, 2000, Codella 
et al., 2004). Although H2O and Class 11 CH3OH masers are produced during 
the early stages of massive star formation, it appears that in general they cannot 
be used to derive the underlying physical processes of the regions. This state- 
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ment does not exclude the possibility that in selected sources very high-spatial 
resolution and proper motion studies of masers can constrain the processes trig- 
gering the maser emission (for prominent recent examples see, e.g., Torrelles 
et al., 2003, Pestalozzi et al., 2004). 

Millimeter continuum emission. Observations of HMPOs in (sub)mm 
continuum emission have shown that they exhibit strong dust continuum emis- 
sion from massive dust and gas cores (e.g., Molinari et al, 2000, Beuther et al, 
2002b, Mueller et al, 2002, Williams et al, 2004). Even with the rather coarse 
spatial resolution of most single-dish investigations (> 11" beam), the gas 
cores often show filamentary and multiple structures. The density distribution 
of the cores are usually well fitted by power-law distributions and the power- 
law indices are in most cases around 1.6 (e.g., Beuther et al., 2002b, Mueller 
et al., 2002), similar to the power-law indices of young low-mass cores (e.g., 
Motte et al., 1998). Furthermore, multi-wavelength investigations show that 
the dust opacity index [3 toward HMPOs is significantly lower than toward 
UCHu regions (mean values of 0.9 versus 2.0, Williams et al, 2004, Hunter, 
1997). In low-mass cores, low dust opacity indices are attributed either to in- 
creasing opacity or grain growth at early evolutionary stages (e.g., Hogerheijde 
and Sandell, 2000, Beckwith et al., 2000). In contrast, Williams et al. (2004) 
show for their sample that the lowest observed values for [3 are associated with 
lower rather than higher opacities. Thus, these observations are consistent with 
grain growth in the early evolutionary stages of the dense cores. However, a 
point of caution remains because significant grain growth is expected to take at 
least 10 5 years (e.g., Ossenkopf and Henning, 1994) whereas the whole mas- 
sive star formation process might only last for ~ 10 5 years (McKee and Tan, 
2002). If these time-scales are correct it may be difficult to produce similarly 
large grains in high-mass cores compared to their low-mass counterparts and 
the observed low values of (3 might require a different interpretation. 

Fragmentation and the Initial Mass Function (IMF). One important ques- 
tion for the whole cluster formation process is: At what evolutionary stage 
does the Initial Mass Function gets established? Do the initial fragmentation 
processes of the massive cores already show signatures of the IMF, or is com- 
petitive accretion during the ongoing cluster formation the main driver of the 
IMF (e.g., Bonnell et al., 2004)? Within the last two years, there have been 
published a few very different studies which investigate the protocluster mass 
distribution and its connection to the IMF: Shirley et al., 2003, analyze single- 
dish CS observations toward a sample of 63 massive H2O maser sources with 
~ 30" resolution, Williams et al., 2004, do a similar study based on single-dish 
(sub)mm continuum observations with 8" — 15" resolution, and Beuther and 
Schilke, 2004, observe one HMPO in the mm continuum emission at 1.5" x 1" 
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Figure 1. Compilation of mass distributions from very young high-mass star-forming regions: the 
cumulative mass distribution to the left is from single-dish CS observations of a sample of massive H2O 
maser sources (Shirley et al., 2003, LS and RE correspond to Least Square and Robust estimation fits). The 
cumulative mass spectrum in the middle is from single-dish mm dust continuum observations of a sample 
of HMPOs (Williams et al., 2004, the dotted and full lines correspond to far and near kinematic distances, 
respectively). The protocluster mass function to the right is from high-spatial resolution PdBI observations 
of one HMPO IRAS 19410+2336 (Beuther and Schilke, 2004, the full line is the best fit to the data, and the 
dashed and dotted lines present the IMFs derived from Salpeter, 1955, and Scalo, 1998). 

resolution with the Plateau de Bure Interferometer. The main difference be- 
tween the two single-dish studies and the interferometric investigation is that 
the single-dish observations average out whole young clusters and thus need 
many clusters to derive a cumulative mass distribution of the sample of proto- 
clusters, whereas Beuther & Schilke (2004) only investigate one young cluster 
but are capable of resolving many sources within this cluster. All of the above 
mentioned studies preferentially select regions that are in the earliest evolution- 
ary stages of star formation and most have not yet formed detectable UCH11 
regions. Figure 1 summarizes the results from the three studies. 

Since the observations and the analysis are very different, the completeness 
limits of the Shirley et al. and Williams et al. studies are different as well (> 
1000 and > 10 M Q , respectively). Nevertheless, both studies find above their 
completeness limits power-law distributions comparable to the IMF. Similarly 
but on much smaller spatial scales, Beuther & Schilke (2004) also find a good 
correspondence between the protocluster mass spectrum in IRAS 19410+2336 
compared with the IMFs as derived by, e.g, Salpeter, 1955, Scalo, 1998. These 
independently derived results are encouraging support for the idea that the 
early fragmentation of massive star-forming cores is responsible for the deter- 
mination of the IMF, and not processes which might take place at evolutionary 
later stages like competitive accretion. An additional result from the single- 
dish observations alone is that the star-formation efficiency for the sample of 
protoclusters is relatively mass invariant (Williams et al., 2004). 

Chemistry. Young massive star-forming regions and especially Hot Cores 
are known to have an incredibly rich chemistry (e.g., the rich spectra toward 
the Orion-KL Hot Core, e.g., Blake et al., 1987, Schilke et al., 1997). Here, 
we present recent submillimeter observations of the prototypical Hot Core 
Orion-KL. A more detailed discussion of core chemistry is presented in Paola 
Casein's contribution to this volume. 

Beuther et al., 2004b, used the Submillimeter Array (SMA) to image the 
submm line and continuum emission of Orion-KL at sub-arcsecond resolution. 
In the submm continuum, the enigmatic source "I" could be disentangled from 
the Hot Core, and source "n" was detected for the first time shortward of 7 mm 
(Beuther et al, 2004b). Additionally, within one spectral setup of the SMA 
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(4 GHz combining the upper and lower sideband) approximately 150 molecu- 
lar lines were observed (Beuther et al. in prep.). Figure 2 shows example line 
images which highlight the spatial variations observed between the different 
molecular species: SiO is detected mainly from the outflow emanating from 
source "I", but all other molecules peak offset from source "I". Typical hot 



Figure 2. Submm line images obtained with the SMA at 865 fim (Beuther et al. in prep.). The grey-scale 
and contours show the molecular emission as labeled at the bottom left of each panel. The three stars mark 
the positions of source "I", the Hot Core mm continuum peak, and source "n" (labeled at the bottom-right). 



core molecules like CH3CN or CH3CH2CN show the hot core morphology 
previously reported from, e.g., NH3 observations (Wilson et al., 2000). Con- 
trary to this, CH3OH does show some hot core emission as well, but there is 
an additional strong CH3OH peak further south associated with the so called 
compact ridge emission. Other molecules like C 34 S or C2H5OH show even 
different spatial emission, tracing partly the hot core and the compact ridge 
but also showing an additional peak to the north-west that is spatially associ- 
ated with IRC6. We still do not properly understand the underlying reasons 
for this spatial molecular diversity, but these observations of the closest Hot 
Core, Orion-KL, stress that we must be cautious in our interpretations of spa- 
tial structures for other Hot Cores that are typically at much greater distances. 

1.2 High-Mass Starless Cores (HMSCs) 

Prior to the infrared satellites ISO and MSX, there was no observational 
evidence for High-Mass Starless Cores, and HMPOs were considered to be 
the youngest observationally accessible stage of massive star formation. The 
detection of Infrared Dark clouds (IRdCs) and their association with mm dust 
continuum emission changed that situation (e.g., Egan et al., 1998), and since 
then observational studies of these HMSCs/IRdCs have begun to shed first light 
on this evolutionary stage. Evans et al., 2002, outlined potential properties of 
HMSC, and recently Garay et al., 2004, presented an analysis of four such 
sources. Here, we present an approach to identify and study HMSCs based on 
a comparison of previously derived mm continuum maps of HMPOs (Beuther 
et al., 2002b) with the MSX database (Sridharan et al. in prep.). 

Based on an initial study of HMPOs associated with IRAS sources (Sridha- 
ran et al., 2002), Beuther et al., 2002b, observed serendipitously within the 
same fields of view a number of secondary mm-peaks which were not as- 
sociated with IRAS sources. Correlation of these mm-peaks with the MSX 
mid-infrared observations revealed a number of mm-peaks which are not only 
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Figure 3. The grey-scale shows the 8 /im emission from the MSX satellite, and the contours present the 
1.2 mm continuum emission first presented in Beuther et al., 2002b. 



weak in the mid-infrared but which are seen as absorption shadows against the 
galactic background (Fig. 3, Sridharan et al. in prep.). 

Observing these HMSCs with the Effelsberg 100 m telescope in NH3 (Srid- 
haran et al. in prep.), the derived average temperatures are around 15 K com- 
pared to the average temperatures of the more evolved HMPOs of 22 K. Fur- 
thermore, we find on average a smaller linewidth of 1.6km/s toward the HM- 
SCs in comparison to a mean value of 2.1 km/s observed toward the HMPOs. 
With gas masses of the order a few hundred M Q , no mid-infrared emission, 
cold temperatures and less turbulence implied by the small line-widths, these 
sources are promising candidates for being HMSCs on the verge of forming 
massive clusters. Ongoing single-dish and interferometric follow-up observa- 
tions of this and other samples will hopefully improve our understanding of 
the initial conditions of massive star formation significantly. 

2. Massive molecular outflows: an evolutionary scenario 

In contrast to low-mass outflows, there is no general census for the under- 
standing of massive molecular outflows. For recent reviews on different as- 
pects of massive molecular outflows see Bachiller and Tafalla, 1999, Church- 
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well, 1999, Richer et al., 2000, Shu et al., 2000, Konigl and Pudritz, 2000, 
Shepherd, 2003, Beuther, 2004. 

Two competing scenarios have been put forth to explain the formation of 
massive stars: one possibility is to extend the low-mass models to high-mass 
star formation and form massive stars via enhanced accretion processes through 
accretion disks (e.g., Jijina and Adams, 1996, Norberg and Maeder, 2000, 
Yorke and Sonnhalter, 2002, McKee and Tan, 2003). Contrary to this, at the 
dense centers of evolving massive clusters, competitive accretion could be- 
come dominant and it might even be possible that (proto)stellar mergers occur 
and form the most massive objects via coagulation processes (e.g., Bonnell 
et al, 1998, Stahler et al, 2000, Bally, 2002, Bonnell et al, 2004). Massive 
molecular outflows can help to differentiate between these scenarios: while the 
enhanced accretion model predicts collimated molecular outflows with proper- 
ties similar to those for low-mass outflows, it is unlikely that highly-collimated 
structures could exist in the coalescence scenario and the outflow energetics 
would be significantly different. 

2.1 Summary of observational constraints 

Observations of massive molecular outflows are unfortunately not conclu- 
sive. On the one hand, some observations find that molecular outflows from 
high-mass star-forming regions appear less collimated than outflows from low- 
mass regions (e.g., Shepherd and Churchwell, 1996a, Ridge and Moore, 2001, 
Wu et al., 2004). However, it has been shown that the lower collimation factors 
of massive outflows observed with single-dish telescopes could simply be due 
to the on average larger distances of high-mass star-forming regions (Beuther 
et al, 2002c). A few interferometric observations with 0.1" — 10" resolution 
suggest that outflow opening angles from at least some young, early-B stars 
can be significantly wider than those seen toward low-mass protostars (e.g., 
Shepherd et al., 1998, Shepherd et al., 2001, Shepherd et al., 2003). Thus, 
greater distances and the associated lower angular resolution is not the only 
reason for the observed poor collimation in some massive outflows. There is 
also evidence that the energetics in at least one massive outflow from an early 
B star in W75N are different to those in low-mass outflows, although exactly 
what this implies is unclear (Shepherd et al., 2003). 

In recent years a few examples of collimated massive outflows from early B 
stars have been found (e.g., Beuther et al., 2002a, Gibb et al, 2003, Su et al., 
2004). A comparison of position-velocity diagrams from collimated massive 
molecular outflows with their low-mass counterparts reveal similar signatures 
from outflows of all masses (Beuther et al, 2004a). Moreover, a near-infrared 
analysis of a HMPO found collimated outflows with properties appearing to be 
scaled-up versions of low-mass outflows (Davis et al., 2004). This latter group 
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of observations favors an accretion scenario for the formation of massive stars 
that is a scaled up version of the low-mass scenario. 

For both well-collimated and poorly collimated outflows, relations between 
outflow and source parameters hold over many orders of magnitude (e.g., Richer 
et al., 2000, Beuther et al., 2002c, Wu et al., 2004). Richer et al. point out that 
these correlations do not necessarily argue for a common entrainment or driv- 
ing mechanism for sources of all luminosities given that one would expect that 
most physically reasonable outflow mechanisms would generate more power- 
ful winds if the source mass and luminosity are increased. It is interesting to 
note that no well-collimated outflows have been observed toward a (proto)-O 
star (e.g., Shepherd, 2003, and references therein, Sollins et al., 2004). Thus, 
the link between accretion and outflow is not as well established for forming 
O stars as for early B stars and formation by coalescence may still be a possi- 
bility. Table 1 presents a number of relevant references for observations which 
shed light on massive outflows. 

The current situation is that different observations find results which appear 
to contradict each other. In order to resolve this conflict, it is necessary to 
reevaluate the observations. In principle, three possibilities exist to solve the 
discrepancies: (a) some observations are wrongly interpreted, (b) there exist 
physically different modes of massive molecular outflows and thus massive 
star formation, or (c) the observed sources are not directly comparable, e.g., 
there might be evolutionary differences between the observed sources which 
could translate into various evolutionary signatures of the massive outflows. 

It is possible that (a) applies to some of the single-dish results, e.g., the 
initial claim of low collimation for massive molecular outflows can be misin- 
terpreted by the low spatial resolution of the observation, the generally larger 
distances to the sources and uncertainties due to the unknown inclination an- 
gles of the outflows. However, the results based on high-spatial-resolution 
observations are more significant and cannot be discarded in this manner. It 
appears that physically different modes and/or evolutionary differences must 
exist to explain all the observations. While we cannot exclude the possibility 
that different modes of star formation and their associated outflows exist, it 
seems unlikely given the link between accretion and outflow that has been es- 
tablished up to early B spectral types. Thus, here we propose an evolutionary 
scenario for massive outflows which is capable of explaining the morphologi- 
cal and energetic results within a consistent picture of massive protostellar evo- 
lution. This evolutionary scenario is based on the idea that stars of all masses 
form via similar accretion-based processes that are proposed for low-mass star 
formation (e.g., McKee and Tan, 2003). 
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Table 1. Summary of outflow results. This table does not claim completeness but gives an 
overview of relevant publications connected to the question of outflow evolution. 



Source 


Results 


Ref 


Low-spatial resolution single-dish studies 


122 UCHns 


ubiquitous phenomena 


1 


10 UCHns 


50% bipolarity, low collimation, M vs Lboi relation 


2 


3 YSOs 


M vs Lboi still holds 


3 


11 YSOs 


low collimation, no Hubble law, M vs Lboi not correlated 


4 


69 HMPOs 


55% bipolarity, accretion-outflow process suggested 


5 


26 HMPOs 


80% bipolarity, no M vs Lboi but M ou t vs M CO re, consistent 
with collimated outflows, similar to low-mass outflows 


6 


139, all types 


M out & M correlate with Lboi, on average less collimated 


7, 8 


High-spatial resolution interferometer and infrared studies 


G 192. 16, HCHll 


wide-angle wind + molecular flow & disk 


9,10,11 


IRAS20126, HCHll 


precessing jet with disk 


12,13 


IRAS05358, HMPO 


consistent with low-mass jets 


14 


W75N, HC/UCHn 


energetics different to low-mass, not scaled up outflows 


15 


IRAS 19410, HMPO 


consistent with low-mass outflows 


16 


G35.2, HCHll 


a cluster of collimated outflows 


17 


2 HMPOs 


accretion based outflows 


18 


4 HMPOs/HCHns 


same physical driving mechanisms as low-mass outflows 


19 


IRAS 16547, HMPO 


collimated jet, similar to low-mass jets 


20,21 


IRAS18151, HMPO 


scaled-up versions of low-mass jets driven by disk accretion 


22 


G5.89, UCHn 


multiple wide-angle outflows, harbors 05 star 


23 



The listed results just refer to mm studies of thermal lines and infrared studies. Ref: (1) Shepherd and 
Churchwell, 1996b, (2) Shepherd and Churchwell, 1996a, (3) Henning et al., 2000, (4) Ridge and Moore, 
2001, (5) Zhang et al., 2001, (6) Beuther et al., 2002c, (7) Wu et al., 2004, (8) Wu et al., 2005, (9) Shepherd 
et al., 1998, (10) Devine et al., 1999, (11) Shepherd et al., 2001, (12) Cesaroni et al., 1999, (13) Shepherd 
et al., 2000, (14) Beuther et al., 2002a, (15) Shepherd et al., 2003, (16) Beuther et al., 2003, (17) Gibb et al., 
2003, (18) Su et al., 2004, (19) Beuther et al., 2004a, (20) Garay et al., 2003, (21) Brooks et al., 2003, (22) 
Davis et al., 2004, (23) Sollins et al., 2004 

2.2 A potential evolutionary scenario 

Reevaluation of the high-spatial resolution literature presented in Table 1 
shows that most of the youngest early B protostars (HMPOs and those produc- 
ing hypercompact Hn regions) support the idea that these high-mass outflows 
are driven by similar physical processes as those seen in low-mass outflows. 
Contrary to this, three sources in particular show different observational char- 
acteristics: G5.89: which harbors an 05 star and multiple wide-angle outflows, 
and G192.16 and W75N: two early B stars with UCHn regions. G192.16 and 
W75N VLA2 both appear to be relatively old (a few x 10 5 years based on the 
CO dynamical timescale), and they have had adequate time to reach the main 
sequence and produce significant ionizing radiation. The outflows are poorly 
collimated within 100 AU of the central star and their outflow morphology on 
larger scales are consistent with wind-blown bubbles. Furthermore, it should 
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be noted that the single-dish outflow studies, which report a lower degree of 
collimation for massive outflows preferentially observed sources with UCHu 
regions (Shepherd and Churchwell, 1996a, Ridge and Moore, 2001). In con- 
trast, most other single-dish surveys were directed toward samples of younger 
sources, namely HMPOs, observing higher collimation degrees consistent with 
low-mass outflows. 

To date, the different observations have been interpreted either in favor of: 
(1) similar driving processes for outflows of all masses (e.g., Beuther et al., 
2002c, Davis et al, 2004); or (2) that massive stars form via different physi- 
cal processes like explosive coalescence or that the outflows are powered by a 
combination of accretion and deflected outflow, (e.g., Churchwell, 1999, Bally, 
2002). Here, we propose a different, less controversial interpretation of the 
data: an evolutionary sequence for massive protostars. We consider two pos- 
sible evolutionary sequences which could result in similar observable outflow 
signatures as shown in Figure 4. 



HMPO 



HCHII 



UC HII 



i \ 



B5-B2 



Bl-08 



Early O 



Figure 4. Sketch of the proposed evolutionary outflow scenario. The three outflow morphologies can 
be caused by two evolutionary sequences: (top) the evolution of an early B-type star from an HMPO via a 
HCHu region to an UCHn region, and (bottom) the evolution of an early O-type star which goes through 
B- and late O-type stages (only approximate labels) before reaching its final mass and stellar luminosity. 

(a) The outflow evolution of an early B-type star: In the accretion scenario, 
a B star forms via accretion through a disk. During its earliest HMPO phase 
no HCHu region has formed yet, and the disk outflow interaction produces 
collimated jet-like outflows (Fig. 4 (left), e.g., IRAS05358). At some point, a 
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HCHu forms and the wind from the central massive star produces an additional 
less collimated outflow component. At that stage the disk is not entirely de- 
stroyed yet, and jet and wind can co-exist (Fig. 4 (center), e.g., IRAS 20126). 
Evolving further, a typical ultracompact Hn region forms above and below the 
accretion disk and the massive wind begins to dominate the whole system (Fig. 
4 (right), e.g., W75N). Remnants of the initial jet might still exist, and density 
gradients within the maternal core are likely to maintain at least a small degree 
of collimation at that stage. 

(b) The outflow evolution of an early O-type star: An intriguing aspect of the 
outflow studies is that so far extremely collimated jet-like outflows have not 
been observed toward sources earlier than Bl. It is possible that this is simply 
a selection effect because the evolutionary timescale of early O stars is ex- 
pected to be only around 10 5 years (McKee and Tan, 2002). Given the scarcity 
of O stars and the short formation time, it is difficult to find sources while they 
are still in the earliest stages of formation. However, to form massive stars 
via accretion, the protostellar objects must accrete even after the central object 
has reached the main sequence. In this scenario, at the beginning, the mas- 
sive protostar reaches a mid-B protostar-like state where it drives collimated 
outflows (Fig. 4 (left)). As the object continues to accrete it becomes an early- 
B to late-0 star, forming a HCHu region while the wind and jet components 
could co-exist (Fig. 4 (center)). Finally, as the central star evolves, it may 
eventually become a mid- to early-0 star. The increased radiation from the 
central star would generate significant Lyman-a photons and ionize the molec- 
ular outflowing gas even at large radii. The result may be an increase in the 
plasma pressure at the base of the outflow which could overwhelm the colli- 
mating effects of a magnetic field (e.g., Shepherd, 2003, Konigl and Pudritz, 
2000). Incorporating the additional physical effects of high-mass protostars 
on the surrounding disk and envelope is one of the challenges of theoretical 
massive outflow modeling (e.g., Konigl, 1999). In this scenario, it would be 
intrinsically impossible to ever observe collimated jet-like outflows from very 
young early O-type (proto)stars. 

2.3 Discussion and potential caveats 

This scenario is not the only possible interpretation, and different physical 
processes could take place at the earliest stages of massive star formation. This 
reevaluation is based only on imaging data of thermal lines in the mm and in- 
frared bands. However, many studies of outflows have also been undertaken by 
the means of maser studies, especially H2O maser. While some maser studies 
found H2O maser emission consistent with jet-like outflows (e.g. Kylafis and 
Pavlakis, 1999), other studies find shell-like expanding features which can- 
not be explained by collimated jets (e.g., Patel et al., 2000, Torrelles et al., 
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2001, Torrelles et al, 2003). The most intriguing case is the H2O maser study 
in W75N: Torrelles et al. (2003) find on projected spatial scales of 1400 AU 
that the two cm continuum sources VLA1 and VLA2 drive very different types 
of outflows, the first one appears to be a collimated jet whereas the second one 
rather resembles a wind-like shell. As the two sources are part of the same 
core, environmental properties cannot explain the difference. Both sources ap- 
pear to have a similar luminosity (B1-B2 spectral types) although the spectral 
type derived for VLA1 should be considered an upper limit due to likely strong 
contamination by ionizing flux produced by shock waves in the jet (Shepherd 
et al., 2004). Assuming the VLA1 jet is produced by an early B star, then the 
difference in morphology would not be a question of the mass of the central 
source as well. Thus, Torrelles et al. (2003) propose that an evolutionary ex- 
planation is most likely. They cannot determine which source is younger but 
based on the extent of the maser emission - the non-collimated features are 
more compact - they suggest that the collimated jet might belong to the more 
evolved source. This would be a counter-example to the scenario presented 
in this paper as well as to the scenarios discussed for low-mass outflows (e.g., 
Andre et al, 2000). However, maser emission appears to be very selective, and 
many known outflow sources have no H2O maser emission at all (e.g., Beuther 
et al., 2002c). Therefore, it is also possible that the collimated part of the out- 
flow from VLA2 is simply not depicted by the data, and the non-collimated 
features have been excited more recently by the central HCH11 region. Fur- 
thermore, the cm continuum emission from VLA1 may be due primarily to the 
jet (Shepherd et al., 2004). Thus, it is feasible that the central source of VLA1 
is either less massive than VLA2 or it is younger and has not formed a HCH11 
region yet. Therefore, we might only observe a collimated feature toward this 
source. The latter scenario is consistent with the evolutionary sequence pre- 
sented here. To properly differentiate between both scenarios in this source 
requires that we find other means to determine the age of the powering sources 
VLA1 and VLA2. 

It should be stressed that the proposed evolutionary scenario is still only a 
potential, qualitative explanation for the observed outflow features which must 
be tested against both theory and observations. 
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